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ABSTRACT 

Convective driving, the mechanism originally proposed by Brickhilj (ll99U [l983) for pulsating white dwarf 
stars, has gained general acceptance as the generic linear instability mechanism in DAV and DBV white dwarfs. 
This physical mechanism naturally leads to a nonlinear formulation, reproducing the observed light curves of 
many pulsating white dwarfs. This numerical model can also provide information on the average depth of a 
star's convection zone and the inclination angle of its pulsation axis. In this paper, we give two sets of results 
of nonlinear light curve fits to data on the DBV GD 358. Our first fit is based on data gathered in 2006 by the 
Whole Earth Telescope (WET); this data set was multiperiodic, containing at least 12 individual modes. Our 
second fit utilizes data obtained in 1996, when GD 358 underwent a dramatic change in excited frequencies 
accompanied by a rapid increase in fractional amplitude; during this event it was essentially monoperiodic. We 
argue that GD 358's convection zone was much thinner in 1996 than in 2006, and we interpret this as a result 
of a short-lived increase in its surface temperature. In addition, we find strong evidence of oblique pulsation 
using two sets of evenly split triplets in the 2006 data. This marks the first time that oblique pulsation has been 
identified in a variable white dwarf star. 

Subject headings: convection — stars: magnetic field — stars: oscillations — stars: variables: general — stars: 
individual (GD 358) 



1. ASTROPHYSICAL CONTEXT 

White dwarf stars offer several advantages for astrophysical 
study. First, they are the evolutionary endpoint of about 97% 
of all stars and are therefore representative of a large fraction 
of the stellar population. Second, the source of their pr es- 
sure support is electron degeneracy ( Chandrasekhar 19391) so 
their bulk mechanical structure is well understood. Third, 
nuclear reactions, if any, contribute a negligible amount to 
their energy, s o their evolution is dominated by simple cooling 
(lMestelll952l) . Finally, they are observed to pulsate in specific 
temperature ranges. The pulsators are believed to be typical in 
every other way, so what we learn asteroseismologically about 
the m should apply to all white dwarf stars (for recent re views, 
see IWinget & Keplerll2008t iFontaine & Brassardll2008l) . 

In addition to learning about the stars themselves, the rel- 
ative simplicity of white dwarfs makes them ideal laborato- 
ries for testing and constraining poorly-understood physical 
processes. One such process, convection, is an important 
energy transfer process in most stars, yet it remains one of 
the largest sources of uncertainty in stellar modeling. For in- 
stance, main-sequence stars at least 20% more massive than 
the Sun have convective cores, and the amount of convec- 
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tive overshoot and mixing is the primary f actor that deter- 
mines their main sequence lifetimes (see, e.g jDi Mauro et aTJ 
120031) . In addition, red giants and AGB stars have large con- 
vective envelopes, and the details of convection play a role in 
the evoluti on of their surface abundances and in their overall 
evolution (iBertelli et al.ll2009h . 

We have developed a method which uses the pulsations 
of white dwarf stars to measure fundamental parameters of 
their convection zones. The physical idea is that the pulsa- 
tions cause local surface temperature ("T e s") variations that 
lead to local variations in the depth of the convection zone. 
As the convection zone waxes and wanes it both absor bs and 
releases energy, modulating t he local energy flux (Brickhill 
119911: iGoldreich & Wulll999l) . Due to the extreme temper- 
ature sensitivity of convection, finite ampl itude pulsations 
can lead to highly nonlinear li ght curves (Bric khiiil [19921 ; 
IWull200 U Using & Koesteiil2001l) . In [Montgomery! (120051) we 
showed how a simple numerical model could be used to ob- 
tain not only good light curve fits but also information on the 
average depth of a star's convection zone and the inclination 
angle of its pulsation axis. 

2. NONLINEAR LIGHT CURVE FITS 



Montgomery (2005) demonstrated that by considering the 
nonlinear response of the convection zone (Brickhill 1992; 
IWull200lh one could obtain excellent fits to the light curves 
of two (nearly) single mode white dwarf pulsators. We 
have since extende d this technique to multi-periodic stars 
dMontgomery|2007l) and have taken into account the nonlinear 
relationship bet ween the b olometric and observed flux varia- 
tions (Montgomery 2008). In this section we describe these 
effects and show how they have been added to our nonlinear 
light curve fitting technique. 

O ur approach i s well- summarized in Montgomery (2005) 
and iMontgomervl (12008). Briefly, we assume that the con- 
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vective turnover time of fluid elements in the convection 
zone is short (< 1 s) compared to the periods of pulsa- 
tion (> 100 s), so that the convection zone responds almost 
instantaneously to the pulsations and is always in hydro- 
static equilibrium. In addition, we assume that the pertur- 
bation of the flux at the base of the convection zone is si- 
nusoidal. Next, we assume that the luminosity changes are 
due only to temperature changes an d not to geometric effects; 
Robinson, Kepler, & Nather ( 1982) found that the fractional 
radius change AR/R ~ 10~ 5 -10~ 4 for the DAVs, with temper- 
ature variations (and the associated changes in limb darken- 
ing) accounting for well over 99 % of the luminos ity variation, 
and this result was confirmed by Watson (1988) for both the 
DAV and DBV stars. Finally, the surface convection zones of 
pulsating white dwarfs are quite thin geometrically, of order 
10~ 5 of the radius of the star; this allows us to use the plane 
parallel approximation and neglect horizontal energy trans- 
port. With these assumptions we are able to derive a rela- 
tionship between the local flux at a given (9, cf>) entering the 
convection zone at its base, fb ase , an d tnat leaving it at the 
photosphere, Fp^: 



Limb-darkening, lo 



phot 



: Fbase + Tc- 



dF, 



phot 



dt 



(1) 



where the new timescale r c = T C (F phot ) describes the changing 
heat capacity of the convection zone as a function of the local 
photospheric flux. We parameterize it as 



TC = T 



T e ff 



'eff.O 



(2) 



where to is the equilibrium value of tc, T e ff is the instanta- 
neous effective temperature and r e ff.o is its equilibrium value, 
and N is a parameter describing the sensitivity of tq to 
changes in T e ff. From standard mixing length theory of con- 
vection we expect that N ^ 90 for DAVs and N ^ 23 for 
DBVs. It is this extreme temperature sensitivity which is re- 
sponsible for the large nonlinearities seen in white dwarf pul- 
sations. For reference, this timescale is closely related to the 
standard thermal timescale (t ±) at the base of the con vection 
zone: for the DAVs r c « At^ dGoldreich & Wull 19991) and for 
the DBVs rc « 0.6 t±. With the further assumption that the 
angular dependence of F\,. ise is given by a spherical harmonic 
Ye m , we can calculate the bolometric flux changes at the sur- 
face of the model, and average them appropriately over the 
visible disk of the model. 

2.1. Improvements in the Modeling 

Since 2005 we have made important technical improve- 
ments to the light curve fitting code. First, we extended the 
code to include the more common multiperiodic case, where 
many modes with different £ and m values are simultaneously 
present. Thus, the flux at the base of the convection zone is 
now given by a sum over the modes: 



Re 



M 

Y^Aje^'^Y^ie^) 

7=1 



(3) 



In this formula, Aj, luj, <j)j, £,-, and ntj are the amplitude, an- 
gular frequency, phase, £, and m values of the j-th mode, and 
the total number of modes is M. 

Second, we adapted the code to simultaneously fit an arbi- 
trary number of observations ("runs"). This is a necessary 
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FIG . 1 . — An example of the limb-darkening in DB model atmospheres with 
logg = 8.0, for the indicated range of temperatures. 



step for applying this technique to multiple runs obtained 
from Whole Earth Telescope (WET) campaigns as well as 
successive nights of single-site observations. Since the code 
allocates and deallocates memory as needed it typically uses 
only 8 MB of RAM, independent of the number of runs in- 
cluded in the fit. 

Third, we replaced our simple analytical prescription for 
limb darkening with tabulated values based on a grid of 
our mod el atmospheres (for a description of the models, see 
lKoesterll20l0h . This grid ranges in r eff from 20,000 K to 
30,000 K in steps of 500 K and in logg from 7.5 to 8.5 in 
steps of 0.25. From this grid we instantaneously calculate the 
local flux as a function of r e ff and /i = cos#; thus, variations 
in the limb darkening with T e ff are automatically included. 
In Figure Q] we show examples of the limb darkening for a 
logg = 8.0 He atmosphere white dwarf model as a function of 
[i for a range of T e fs values. 

Finally, we improved the way in which the local bolometric 
flux variations are mapped into variations in a given wave- 
length band. Previously we used a flux "correction factor," 
ax, to accomplish this. Denoting by Fx the flux in the pass- 
band X, then, for small fractional changes in the fluxes, ax 
was defined by 



SFx 
Fx 



5E 



■ a x - 



phot 



"phot 



(4) 



where f p h ot is the local bolometric flux at the photosphere and 
SF is the variation in the respective fluxes due to the pulsa- 
tions. Clearly, ax depends on the wavelength coverage of 
the passband as well as the wavelength distribution of the 
flux from the source. In previous analyses we estimated that 
a ~ 0.42 for DBVs and a ~ 0.66 for DAVs. The value for 
the DAVs is not that different from what one obtains from a 
proper calculation assuming a passband centered on 5000 A. 
For the DBVs, however, the more detailed calculations yield 
a value of a ^0.25-0.35 which is significantly different from 
the earlier estimates. 
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FIG. 2. — Top panel: the bolometric flux (F p hot) and the flux in the given 
passband (Fx ) as a function of 7 c ff for fj, = 1 . These calculations are for a 
logg = 8.0 DB model, and both fluxes have been normalized to one at T e g = 
25,000 K. Lower panel: the passband flux as a function of the bolometric 
flux, again normalized at T t ({ = 25,000 K. This is the function / given by 
equation \5\ 

Of course, a is not strictly a constant but rather is a function 
of J e ff and therefore F p h t- For larger amplitudes, departures 
from linearity between the fluxes become more important, and 
to do the problem properly we need Fx as a function of Fp no t, 
i.e., F x = Fx(Fphot)- 

To calculate this flux conversion we use the same grid of DB 
model atmospheres described earlier. In each model atmo- 
sphere, the flux is tabulated as a function of wavelength A and 
viewing angle \i, so limb darkening is automatically included. 
Using a given passband, we integrate it against the model at- 
mosphere flux to obtain Fx as a function of T & s, fi, and \ogg. 
Since logg is essentially constant during the pulsations, we 
only have to interpolate once between sets of models for the 
given value of logg; thus, we have Fx = Fx(T e fs,/j,). Finally, 
since F p hot oc T* s , we can express this as Fx = Fr(F p h t, 

In practice, we are only interested in relative flux changes, 
so we calculate everything relative to a reference flux, which 
we take to be the equilibrium flux of the star. If we denote by 
a subscript "0" the equilibrium values of the fluxes, then the 
function / we want is defined by 

Fx ,( .Fphot \ ... 

p = ■ \ ~f — >H- (5) 

fX,0 V^phot.O / 

In the top panel of Figure [2] we show the fluxes Fx /Fx ,o 
and F p hot /^phot.o as a function of r ett ', both normalized to one 



TABLE 1 

Observing Runs Used for Light Curve Fits 



Run Name 


Telescope 


Instrument 


Length 


chin20060527 


BAO 2.16m 


PMT 


4.5 


chin20060528 


BAO 2.16m 


PMT 


5.2 


chin20060531 


BAO 2.16m 


PMT 


3.9 


hawa20060518 


0.6m 


Apogee 


2.2 


hawa20060519 


0.6m 


Apogee 


1.7 


hawa20060520 


0.6m 


Apogee 


3.5 


hawa20()60521 


0.6m 


Apogee 


7.8 


hawa20060522 


0.6m 


Apogee 


6.6 


hawa20060523 


0.6m 


Apogee 


3.7 


hawa20060524 


0.6m 


Apogee 


5.7 


hawa20060525 


0.6m 


Apogee 


8.8 


hawa20060526 


0.6m 


Apogee 


9.1 


hawa20060527 


0.6m 


Apogee 


7.8 


hawa20060528 


0.6m 


Apogee 


9.1 


hawa20060530 


0.6m 


Apogee 


8.7 


kpno20060518 


KPN0 2.1m 


Apogee 


7.0 


kpno20060519 


KPN0 2.1m 


Apogee 


7.3 


kpno20060520 


KPN0 2.1m 


Apogee 


7.6 


kpno20060521 


KPN0 2.1m 


Apogee 


7.3 


kpno20060522 


KPN0 2.1m 


Apogee 


1.0 


kpno20060523 


KPN0 2.1m 


Apogee 


2.9 


mcdo20060523 


2.1m 


Argos 


7.4 


mcdo20060524 


2.1m 


Argos 


7.2 


mcdo20060525 


2.1m 


Argos 


6.4 


mcdo20060528b 


2.1m 


Argos 


7.2 


mcdo20060529 


2.1m 


Argos 


8.2 


nord20060607 


2.7m 


ALFOSC 


7.1 


nord20060608 


2.7m 


ALFOSC 


8.0 


nord20060609 


2.7m 


ALFOSC 


7.9 



at r e ff = 25,000 K; the filter X is assumed to have an ef- 
fective wavelength of ~ 5200 A. The response of this filter 
was obtained from convolving the wavelength response of the 
ALFOSC-FASU CCD, the NOT primary and secondary mir- 
ror reflectivities, an S8612 filter, and atmospheric absorption. 
In the lower panel we show the function / defined in equa- 
tion [5] as derived from these calculations. We see that while 
this is not a perfectly linear relationship, the deviations from 
linearity are not dramatic. Thus, while we use the fully non- 
linear relation in our calculations, we expect this effect to 
make only a minor contribution to the overall nonlinearities 
associated with the pulsations. 

3. LIGHT CURVE FITS TO THE 2006 WET RUN 

Our recent work with light curve fitting has been limited 
to nearly single mode pulsators: the DAVs G29-360 and 
GD154, and the DBV PG1351+489. This is because (1) 
mono-periodic data can be folded at the pulsation period, pro- 
ducing a high S/N "light curve," and (2) the number of possi- 
ble mode identifications {£ and m values) for a single mode is 
small enough that all possibilities can be directly explored. 

GD 358 violates both of these conditions. First, due to the 
nonlinear interaction of its large amplitude modes, the pulse 
shape obtained by folding its light curve at a mode period 
is not the same as the pulse sh ape which would b e obtained 
in the absence of other modes (Montgomery 2007). Second, 
GD 358 has a large number of observed modes (see Table [2J 
and it is impractical to search all possible combinations of I 
and m values which each mode can take. For instance, if we 
assume GD 358 to have of order ~ 10 modes, all of which 
have I = 1, then all possible permutations of m values yield 

8 Mor e precisely, while G 29-38 is normally multiperiodic, the data set 
used by Monteomery (2005) was taken at a time when its light curve was 
dominated by a single large mode. 
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FIG. 3. — A comparison of the simultaneous fit of Tablef5](solid line) to the light curve from run mcdo20060523 (crosses). 



a number of (21+ l) 10 ~ 60,000 different cases! Since each 
fit takes of order an hour on a single processor computer, this 
is completely impractical using a standard desktop computing 
approach. 

Fortunately, GD 358 is well studied, so we have a good 
idea what the £ and m values for the main pulsation modes ar e 
(e.g.. iMetcalfe. Nather. & Wrn get 2000: IWinget^tai][T994T) . 
Even so, our derived values of to depend only weakly on the 
assumed mode identifications, so it is not necessary for the 
mode identification of each mode to be exact. Furthermore, 
the extended time baseline and excellent coverage of the 2006 
WET run all ows us to obtain very accurate frequencies for 
these modes dProvencal et al.ll2009l) . With this as a basis, our 
approach is to assume that the frequencies and mode identi- 
fications are known and then to make nonlinear light curve 
fits to a subset of runs in the WET campaign which have high 
S/N. This implicitly assumes that the pulsations are coherent 
throughout the time spanned by the runs; in section|6]we show 
that this is not strictly true for some of the modes. High S/N 
data are desirable since we are mainly interested in the non- 
linear part of the light curve, which itself is smaller than the 
linear part. 

Many high quality runs were taken during the May 2008 
WET run, and we list in Table [T] those used in our fits. We 
included the 12 largest amplitude periodicities from the 2006 
WET run which were deemed to be independent frequencies 
and not linear combinations, and these are given in Table [2] 



TABLE 2 
Independent Modes from 
the 2006 WET Run 



Period 

(s) 


k I 


m 


422.56 


8 


1 


423.90 


8 1 


-1 


463.38 


9 


1 


464.21 


9 





465.03 


9 1 


-1 


571.74 


12 


1 


574.16 


12 1 





575.93 


12 


-1 


699.68 


15 





810.29 


18 





852.50 


19 





962.38 


22 






The £ and m identifications are taken from previous analyses 
of this star dWinget et alJlT994t iKepler et alj|200lt iMetcalfd 
120031: iProvenca let al]l2009ft and are based on asymptotic the- 
ory as well as genetic algorithm fits. To calculate the conver- 
sion from bolometric to observed passband fluxes, we must 
assume values for T e g and logg. We chose the values of 
iBeauchamp et all d 1999ft . r eff = 249 00 K and logg = 7 . 91, al- 
though we have also used those of Castanh eira et al.l (12005b 
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FIG. 4. — The variation of £>,, N, and to as a function of time during the 2006 WET run. While to and N show no obvious signatures, 0, shows variations which 
may be sinusoidal in origin. 



(T e ff = 24100 K, \ogg = 7.91) to assess the uncertainties this 
introduces. 

We began the fitting process with the highest S/N data taken 
with the 2.7m Nordic Telescope (NOT). However, it became 
apparent that GD 358's complex pulsation spectrum required 
an extended baseline of data to constrain the phases of closely 
spaced frequencies. Fortunately, GD 358's brightness and 
large amplitude meant that a large number of individual ob- 
serving runs met the required S/N. 

Our simultaneous fit to the runs in Table [T] yields the fol- 
lowing parameters: ro = 572.9 ±6.1 s, N = 23.5 ±0.1, 0; = 
50.5 ± 0.2°. Figure [3] shows the ability of the fit to reproduce 
the essential features of the light curve (run mcdo20060523 
plotted). Additional results for this fit, including the ampli- 
tude and phase for each mode, are given in Table[3] The given 
error bars are formal and should be treated as lower bounds. 

To test the sensitivity of the value of tq to the £ identifi- 
cations we re-computed fits changing the assumed £ of the 
largest amplitude mode (P = 810.29 s) from £ = 1 to £ = 2. 
With this assumption, the best fit resulted for an identifica- 
tion with I = 2, m = 1 and had the following parameter values: 
r = 569.5 ±7.6 s,N= 14.5 ± 0.1, 6 t = 61.3 ± 0.1°. This indi- 
cates that the value obtained for tq does not crucially depend 
on the mode identifications of each mode in the fit. 

Examination of the fit for each individual run revealed an 
apparent modulation of amplitudes, i.e., on some nights the 
variations in the light curve were smaller than the fit and on 



TABLE 3 

Simultaneous Fit to GD 358 Data Set: 
tq — 572.9 i 6. 1 s, W = 23.5±0.1, = 5O.5±O.2° 



Period t 


m 


Amplitude 


Phase (rad) 


962.385 





0.1087 ±0.0012 


2.4641 ± 0.0069 


852.502 1 





0.1198 ±0.0015 


3.3007 ± 0.0075 


810.291 1 





0.4581 ±0.0049 


3.1301 ± 0.0030 


575.933 1 


-1 


0.4838 ±0.0051 


3.4434 ± 0.0039 


574.162 1 





0.2257 ± 0.0024 


5.5258 ± 0.0062 


573.485 1 





0.1082 ±0.0016 


3.6576 ± 0.0115 


571.735 1 


1 


0.3728 ± 0.0040 


2.4489 ± 0.0046 


465.034 1 


-1 


0.1408 ±0.0021 


1.8310 ±0.0123 


464.209 1 





0.1391 ±0.0018 


2.4672 ± 0.0097 


463.376 1 


1 


0.2540 ± 0.0030 


0.2052 ± 0.0073 


423.898 1 


-1 


0.2406 ± 0.0030 


2.0455 ± 0.0083 


422.561 1 


1 


0.2537 ±0.0031 


0.6020 ± 0.0077 



others they were larger. In addition, looking at data from 
single sites suggested that the change from a smaller to a 
larger amplitude state alternated on a roughly night to night 
timescale. 

To test whether a geometric effect or a change in the back- 
ground state of the star could be causing the night-to-night 
modulation of the amplitudes, we went back and re-fit each 
run in Table [T] individually. For these fits we fixed the pe- 
riod and amplitude of each mode to be those given in Table[3] 
but we allowed the inclination angle 6>, and the phases of the 
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FIG. 5. — A section of the light curve of GD 358 during the sforzando event 
in 1996. The crosses are the data points and the curve is the best nonlinear 
light curve fit. 



modes to vary. We also allowed to and N to vary. Thus, from 
each run we obtained a best fit value of the parameters at the 
time of the run. In Figure [4] we present the results of this pro- 
cedure; we plot the variations in 9,, N, and to as a function 
of time. While no clear trends are seen in To and N (lower 
panels), the variations in 9, are suggestive of a periodic ori- 
gin (top panel). Assuming a sinusoidal variation, we have 
included the best fit sine curve in this plot; it has an ampli- 
tude of 14.9 ± 1 .8° and a period of 2.17 ± 0.01 days. On the 
other hand, if we ignore the possible origin of these variations 
and simply consider them to be separate measures of to, N, 
and 9,, then we obtain fairly conservative limits on the error 
of the mean values of these quantities: to = 586.0 ± 11.7 s, 
N = 21 A ± 0.7, and 0* = 47.5 ± 2.2°. We return to the ques- 
tion of this modulation in section [6] 

4. THE SFORZANDO 

During the 2006 WET run, GD 358 was stable in that the 
amplitudes and phases of its modes were fairly constant over 
the 3 week length of the run (except for the modulation de- 
tailed above and in section [6j. However, GD 358 is known 
to change its pulsation spectrum on a range of timescales. A 
spectacular example of this behavior occurred in 1996. Within 
a period of 36 hrs, all of the power in the high k range (pulsa- 
tions with periods > 700 s) disappeared within detection lim- 
its. At the same time, GD 358 more than doubled its apparent 
pulsation amplitude, with power appearing almost exclusively 
at lower k, with a period of ~ 420 s. Over the next week, 
its amplitude decreased to "normal" levels, while the high 
k power did not return for approximately one month. This 
dramatic change is docume nted in several papers, mo st re- 
cently in lKepleret all (120031) and lProvencal et alJ (120091) . The 
episode itself is termed the "sforzando" after a musical term 
for a sudden and short-lived increase in loudness. 

In Figure|5]we show a section of GD 358's light curve dur- 
ing the sforzando. Ironically, during this dramatic episode 
the light curve was much simpler. While the amplitude was 
larger, the star was essentially a single-mode pulsator and the 
individual pulses were more sinusoidal (less nonlinear) than 
before. Still, GD 358 was in a transient state and not all pulses 
were identical. The upper and lower panels of Figure|5]show 
pulses significantly below and above the fit, while the mid- 
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die panel seems to represent a repeating, "quasi-static" state. 
Even with these caveats, the relative simplicity of these data 
are ideal for nonlinear light curve fitting. 

From a time series analysis of this short section of data, we 
find that these pulses have a period of approximately 420.7 s. 
Assuming £ = 1 and trying all values of to, we find the best 
fit shown in Figure |5] to = 0, t = 41.6 ±2.3 s, N = 3.6 ±0.2, 
and 9j = 56.1 ±1.1. For I = 1, to = 1 a similar quality fit can 
be found, but it has a less plausible inclination angle (~ 85°) 
that is not consistent with the values of 9, previously found in 
this paper; even so, to = 28.1 s for this fit. If we assume 1 = 2 
then the best fit has to ~ 24.2 s and to = 0, although it requires 
such a large intrinsic amplitude that we consider it unphysical. 
Summing up, while we have good reasons for preferring the 
I = 1, to = identification for our fits, the overall value of to 
is not strongly dependent on this identification. 

Comparing these results to those of section [5J we find that 
To was much smaller during the sforzando than it was during 
the 2006 WET run. This implies that, for whatever reason, 
GD 358's convection zone was thinner during the sforzando. 
As we can see from Figure [6j this would imply a r e ff sev- 
eral thousand degrees hotter th an its normal temperatur e. For 
ML2/a = 1 .0 convection (e.g., Bo hm & Cassinelll 19711) . this 
would mean T & s ~ 27,000 K, which is two to three thousand 
degrees hotter than its normal equilibrium state. 

Such an increase in temperature would lead to about a 40% 
increase in the bolometric luminosity. From Figure [2] we see 
that such an increase in the bolometric luminosity translates 
into a 15-20% increase in intensity as measured in a pass- 
band centered at 5200 A. While such a passband is reason- 
able for a CCD + BG40/S8612 filter + atmosphere, the ob- 
servations in 1996 were made with phototubes, and these are 
much more blue sensitive. Thus, for a phototube passband 
centered around 3800 A we expect intensity increases of about 
20-25%. 

As reported in Table 5 of iProvencal et alJ (120091) . such an 
increase in intensity was observed in GD 358 during the 
sforzando event. In Figure|7]we plot these data, where the tri- 
angles are data from McDonald Observatory and the squares 
are data from Mt. Suhora Observatory in Poland. Both curves 
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are normalized to a value of one near BJED ~ 24503 1 1 . We 
see that these data, although sparse, strongly suggest that 
GD 358 underwent a sharp intensity increase and decline of 
about 20-30% during this event. This corroborates our ear- 
lier interpretation that GD 358 had a thinner convection zone 
with a smaller value of to, and that this thinness was due to its 
surface layers being temporarily hotter. 

We note that simil ar results and conc l usions have previ- 
ously been made by iWeidner & Koesterl (120031) . Using the 
approach of Ising & Koester (2001J) they were the first to 
make numerical simulations of GD 358's light curve during 
the sforzando event. As in our present analysis, they found 
that the shape of the light curve could only be reproduced for 
higher effective temperatures than normally assumed for this 
star, suggesting a value of T e g ~ 27,000 K. 

4. 1 . The Amplitude of the k=8 Mode 

The physical origin of the sforzando event is completely un- 
known. All we know for certain is that the amplitudes of the 
high-fe modes rapidly decreased, leaving power in the k = 8 
mode, which itself increased in amplitude by approximately 
a factor of 30. While we have no model to explain the dis- 
appearance of power in the high-A: modes, in this section we 
consider what would happen to the apparent amplitude of the 
k = 8 mode if the convection zone were suddenly removed. 

One effect of the thermal response of the convection zone 
is to attenuate the amplitude of the flux variations incident on 
its base. It acts as a low-pass filter, reducing the photospheric 
amplitude of a mode according to 

AFphot _ 1 Ai*b aS e 

Fa ^/ \ + (utq) 1 Fq 

where Fq is the equilibrium value of the flux, Fb ase and F p h ot 
are the instantaneous fluxes at the base of the convection zone 
and at the photosphe re, respectively, and ui i s the angular fre- 
quenc y of the mode ( Goldreic h & Wulll999t IWu & Gol dreich 
1999). 

If we assume that the amplitude of the k = 8 mode at 
the base of the convection zone was constant throughout the 
sforzando, equation [6] shows that a decrease in to would nat- 
urally lead to an increase in observed amplitude. Thus, for 



a 420 s mode we would see an increase of approximately a 
factor of 8 in its apparent amplitude as to goes from ~ 600 s 
to ~ 40 s. This is a large factor, but it is still much less than 
the factor of ~ 30 increase that was actually observed. Thus, 
the k = 8 appears to have increased its intrinsic amplitude by 
of order a factor of 4 during the sforzando. How it was able 
to do this on so short a time scale remains a mystery. 

5. THE CONVECTIVE TIME SCALE t 

In Figure[6]we plot the known determinations of the convec- 
tive response timescale, to, versus T e s for the DBVs. Based 
on the discussion in the previous section, GD 358 is plotted 
twice: once for the 2006 data and once for the 1996 data. 
In addition, we show the positi on of the DBV PG 135 1+489 
(tq ~ 100 s, Montgomery 2005), where we have assumed the 
pure He fit for its T e g (Beauchamp et al. 1999). As expected, 
the data indicate an increase in the depth/mass of the convec- 
tion zone with decreasing T e s. 

We also plot line s in F igure [6] s howing th e predictions of 
ML2/a convection (Bohm & Cassinelli 1 9711) for various val- 
ues of a. We see that while low values of a are excluded 
(a < 0.6), values in the range 1.0 to 1.2 provide a reasonable 
description of how to varies with 7^. We note that the \ogg 
values determined for these stars are nearly identical, so they 
do form an actual sequence in T e ff. In general, however, to is 
also a function of \ogg, albeit a somewhat weaker one. 

Our ultimate goal is to map to as a function of both T & s and 
logg for both the DBV and DAV instability strips. This will 
provide important reference points for new hydrodynamic 
simul ations of convection which are starting to come online 
(e.g. iMuthsam et aT] 120091) . For instance, for a given white 
dwarf one can use the measured r en - and logg values and 
perform a hydrodynamic simulati on o f its c onvection zone. 
Then, using equations 3 and 5 of iWul ( 120011) one can com- 
pute to. This value of to can then be compared with the value 
derived from the light curve fits outlined above. 

6. THE OBLIQUE PULSATOR MODEL 

6.1. The Formalism 

In section [3] we found evidence of a modulation of the in- 
clination angle 6>, (see Figure @J. This modulation has a for- 
mal significance level of 14.9/1.8 ~ 8<7, which cries out for a 
physical interpretation. The most obvious is some form of the 
"oblique pulsator" model, in which the pulsation axis is in- 
clined with respec t to the rotation axis, and this axis prece sses 
as the star rotates (Kurtz & Shibahashi|[l986t lKurtzlll982l) . A 
magnetic field that is inclined to the rotation axis is usually 
invoked, and the pulsations are assumed to be aligned with 
the magnetic axis. 

While this hypothesis does introduce several unknowns (a 
magnetic field, the angle between the magnetic and rotation 
axes, etc.) it also makes three testable predictions. If we take 
/ and / rot to be the mode frequency in the frame of the star 
and the rotational frequency, respectively, then an £ = 1 mode 
with frequency / aligned with the magnetic axis will appear 
as three separate peaks in the Fourier transform of the light 
curve, with frequencies of /— / rot , /, and f + f - t- For clarity 
we will refer to these peaks as "geometric" peaks as they only 
appear because the pulsation axis spins around the rotation 
axis, leading to a periodic apparent amplitude modulation of 
the mode. This amplitude modulation manifests itself in the 
Fourier transform as two additional "geometric" peaks on ei- 
ther side of the original frequency, with the beating of these 
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TABLE 4 

Frequency solution for oblique pulsation model: 
/ rot = 5.362 ±0.003 fiHz 



Frequency (fiHz) 


Amplitude (mma) 


Phase (rad) 


A$/2tt 


Independent Freqs." 


Triplet 1 


1736.302 ± 0.004 
1741.664 ±0.003 
1747.027 ± 0.004 


16.77 ±0.13 
10.81 ±0.13 
1.75 ±0.13 


0.739 ± 0.008 
2.832 ± 0.012 
1.724 ±0.075 


0.510 ±0.013 


1736.302 ± 0.001 
1741.665 ±0.001 
1746.673 ± 0.007 


Triplet 2 


1738.362 ± 0.005 
1743.725 ± 0.003 
1749.087 ± 0.005 


0.95 ±0.13 
5.56 ±0.13 
12.55 ±0.13 


1.268 ±0.138 
0.174 ±0.023 
2.117 ±0.011 


0.517 ±0.023 


1737.962 ± 0.007 
1743.738 ± 0.002 
1749.083 ± 0.001 



1 These are the unconstrained frequency fits of Provencal et al. 12009). 



three peaks producing the periodic amplitude changes. As this 
splitting is caused by the rotation of the star, these frequencies 
must be equally split to within measurement errors. This is the 
first and most stringent condition the model must face. 

Second, the oblique pulsator model predicts a specific phase 
relationship for the peaks in a given geometric triplet (or 
for hig her £, a (21+ l)-multiplet). Frorn lKurtz & Shibahashi 
(1986), an £ = 1 mode generically has luminosity variations 
given by 



AL /L=A-Cos[(aj-tt)t + 4>] 

+ A cos [tut + (/)] + A+ cos [(w + £l)t + (/)] , 



(7) 



where oj and are 2ir times / and / lot , respectively (see Ap- 
pendix [A] for the complete expressions). These components 
will only have the same phase <p for a particular choice of the 
zero point of time. For other zero points one can show that 
this phase relation translates to 



2$ -($ + +$-) = 0, 



(8) 



where {$_, $o, are the measured phases of the respective 
components. In general, the product A_ x A+ can be negative. 
For this case, if we define the amplitudes always to be positive 
and absorb any minus signs into the phase for each peak, the 
relation becomes 



2$ -(<£++$_) = 7T. 



(9) 



Defining A$ = 2<P - ($+ + $-), then from the sign of the 
amplitudes in equations (IA2b — (IA4t we see that A$ /2tt = 
for m = modes and A<I>/27r = 0.5 for m = ±1 modes. These 
phase relations, while less iron-clad than the equal spacing of 
the triplets, should be satisfied within the errors for geometric 
peaks split by oblique rotation. 

A third condition/prediction is the relative amplitudes of the 
geometric peaks. Such a calculation makes assumptions con- 
cerning the nature and strength of the magnetic field, so this 
prediction of t he mo del is the least reliable of the three. In 
equations ( IA2b — ( IA4b in Appendix lAl we give expressions for 
the amplitudes of the various components of I = 1 modes per- 
turbed by a magnetic field and oblique rotation. The relevant 
parameters are the inclination angle of the rotation axis, 
the obliquity of the magnetic axis, /3, and x\ . The parameter 
x\ is given by 



(10) 



where Cki is the rotational splitting coefficient due to the Cori- 
ohs force, ft is the angular frequency of rotation, and Uq 



and oj' 1 ' are the perturbations to the frequencies of the m = 
and m = 1 intrinsic modes due to the magnetic field, respec- 
tively. Even though this is clearly a more model dependent 
statement, we would still hope that the amplitudes could ap- 
proximately be fit and/or predicted within the oblique pulsator 
formalism. 



6.2. The Results 

First, we consi der a fit to the k= \2 re gion of the FT, from 
1730-1750 /zHz. iProvencal et ail d2009l) found 6 significant 
peaks in this region. We interpret these peaks as resulting 
from 2 components of an £ = 1 mode, each split into a geo- 
metric triplet by oblique pulsation. These two original peaks 
are part of an intrinsic triplet produced by standard rotational 
splitting, but the amplitude of the third member is below our 
detection threshold. 

We fit two sets of exactly evenly split triplets to the data set, 
where the central frequencies of each triplet are 1741.664 and 
1743.725 /iHz, and the value of the splitting is 5.362 fiHz. 
The values of the fit parameters we obtained are given in Ta- 
ble [4] In Figure [8] we show the result of pre- whitening the 
k = 12 region by this solution. The reduction in power of the 
FT is very significant, showing that evenly split triplets are a 
good representation of the data. This is a necessary condition 
for the oblique pulsator model to be applicable. 

Second, we consider the phase relations within each of 
these equally split triplets. From Table|4]we see that the first 
triplet has A$/27r = 0.510±0.013 and that the second triplet 
has A$/27r = 0.517 ±0.023. The oblique pulsator model also 
passes this test with flying colors. In addition, the fact that 
A$/27r ~ 0.5 rather than 0.0 implies that the modes are not 
axi-symmetric, i.e., \m\ = 1 for both modes. For whatever rea- 
son, the m = member of the original triplet is not present at 
observable amplitudes. 

Finally, we wish to test the ability of the oblique pulsator 
model to adequately reproduce the amplitudes of the geomet- 
ric peaks. To calculate the amplitudes w e use the analyti- 
cal expressions given in A ppendix lAl (see lUnno et al] 119891 : 
Kurtz & Shibahashi 1986|) and to perform the fits we have 
used a genetic algorithm dCharbonneaulll99"5l) . This allows 
us to search the m values of each intrinsic triplet as well as the 
values of the parameters (3, and x\ . As we demonstrate be- 
low, such a fit is more constrained than one would think given 
the six data points and five free parameters. 

In Figure [9] we show the result of the fit to the amplitudes in 
Table|4] The fit is quite impressive. It has the added bonus that 
triplets 1 and 2 are required to originate from \m\ = 1 intrinsic 
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FIG. 8. — The effect of pre- whitening by two identically split triplets in the 
oblique pulsator model. The top panel is the original FT of the k=\2 region 
and the lower panel is the result of pre-whitening by the two pairs of triplets. 
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FIG. 10. — Pre-whitening of the same region as in Figure [8] but using the 
full oblique pulsator model for the amplitudes, frequencies, and phases. Thus, 
not only are the frequency splittings constrained, but the amplitudes within a 
triplet and their relative phases are as well. 



modes and they must have opposite signs. This corroborates 
our earlier result based on the phases that both triplets origi- 
nated from \m\ = 1 intrinsic modes. 

Given the near equality of the number of data points and 
free parameters we wished to assess how easily our model 
could reproduce any data set. To do this, we randomly gen- 
erated amplitudes for 500 pairs of triplets and fit them in the 
same way we fit the data, normalizing the residuals by the 
mean squared amplitudes of the two triplets. We found that in 
only 1 1 cases out of 500 were the random amplitudes better fit 
by our model than the measured amplitudes were. Thus, we 
conclude that our amplitude fits are significant at the ^98% 
level. 

As a final check on this procedure, we use these amplitudes 
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amplitudes (points) for the different components of "k=12" assuming the 
oblique pulsator model. The errors on the observed amplitudes are smaller 
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together with equations IA2I and IA4I to pre- whiten the k = 12 
region of the FT, as shown in Figure [10] Although consid- 
erably more constrained, this pre-whitening is clearly just as 
effective as that shown in Figure|8] Thus, the oblique pulsator 
model passes all tests with flying colors when applied to the 
k = 12 region of the FT. 

6.3. Interpretation of the Oblique Pulsator Fits 

The relative amplitudes of the peaks within a triplet split 
by obli que pulsation/rotation depend on f3, and x\ (see 
section 16.11 for a definition of these quantities). From the 
amplitude fits obtained in the previous section, we find that 
6i = 46.3°, p = 31.8°, and x x = 5.65. This value of <?,• is close 
to that obtained from our nonlinear light curve fits, 47.5 ±2.2° 
(section[3]i. This provides an important consistency check on 
both methods. 

Given that is the angle between the rotation axis and our 
line-of-sight and that f3 is angle between the rotation and pul- 
sation/magnetic axes, the effective inclination angle of a mode 
aligned with the magnetic axis should vary between - (3 and 
Oj + P as the star rotates, i.e., between 14.5° and 78.1°. From 
the top panel of Figure|4]we see that 0, varies with about half 
of this amplitude. The straightforward explanation for this is 
that only the k = 12 modes are affected by oblique rotation, 
while the other modes are aligned with the rotation axis. As 
a result, the value of 9, shown in Figure|4]represents the aver- 
age inclination angle for all modes. Averaging over all modes 
yields an inclination angle which varies by about 15° rather 
than the full 32° experienced by the k = 12 modes. 

Finally, the parameter x\ is a measure of the relative 
strength of the Coriolis and magnetic perturbations; the fit 
value x\ =5.65 indicates that the magnetic field is the dom- 
inant perturber. This reinforces the idea that the pulsations of 
these modes are tied to the magnetic axis. In addition, for a 
large-scale dipole field simple perturbation theory applied to 
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' = 1 and 2 modes yields 



(i) 

M i 



| + |m 2 , l=\ 



f-fm 2 , 1 = 2 



(11) 



where wffi is the perturbation to the frequency of the mode 
due t o the magnetic field (Montgomery 1994; [Jones et aTl 
1989). Together with equation ([Toi l, the fact that x\ is posi- 



tive means that w^ 1 ' > uj\ l> , which in turn suggests that £ > 2. 
In reality, this indicates that these modes, at least in the sur- 
face layers, may not be pure £ = 1 modes but contain a mixture 
of other higher £ components. 

The success of the oblique pulsator model in explaining the 
k=\2 peaks constitutes evidence that oblique pulsation plays 
an important role in GD 358. However, there is little con- 
clusive evidence that other modes in this star are undergo- 
ing oblique pulsation. On the contrary, the stable k = 8 and 
9 triplets are consistent only with the traditional case of the 
pulsation axis aligned with the rotation axis. A possible ex- 
planation is that the outer turning points for the k = 8 and 9 
modes are farther from the surface than they are for the k=l2 
modes, so if the dominant magnetic effects are also confined 
to the surface layers then the lower k modes would not be af- 
fected by the magnetic field and the higher k modes would. 

The magnitudes of the multiplet splittings/fine structure 
also present a challenge. For instance, the oblique pulsator 
model requires a rotation frequency of 5.36 /iHz, whereas 
the splittings of the k = 8 and 9 modes, believed to result 
from standard rotational splitting, give a value of 7.5 /iHz 
dProvencal et alj I2009t IWinget et : all 119941) . Since different 
modes sample different regions of the star's interior, differen- 
tial rotation in the radial direction could explain these results. 

A further concern is the internal consistency of the oblique 
pulsator fit. The central components of each of the triplets in 
Table [4] should themselves be members of an intrinsic triplet 
which is split by rotation. This splitting is observed to be 
2.06 /iHz and should equal 2Q£/ rot cos/3 for solid body ro- 
tation of the star. This is satisfied if Cut ~ 0.225 instead of 
Ckt ~ 0.5 as is expected for £ = 1 modes. The splitting within 
a geometric triplet is a direct measure of the surface rotation 
rate, while the difference in frequency of the central com- 
ponents of these triplets (the "intrinsic" modes) is given by 
2C J t^/ rot cos (3, where /, ot represents a radial average of the 
rotation profile. Thus, differential rotation, with the interior 
rotating more slowly than the surface layers, could explain 
these values. An alternate explanation is that the magnetic 
field mixes higher £ components into the spatial structure of 
the k = 12 modes. They would therefore have rotational con- 
stants Cue indicative of higher £ modes. Since Cut ~ !/£(£+ 1) 
for moderate to high k modes, these modes would have corre- 
spondingly smaller values of Cug- 

A further possibility is that the intrinsic modes have £ = 2 
instead of I = 1 . In this case only 3 of the 5 possible geo- 
metric peaks would be large enough to be detectable. First, 
we note that the predicted phase relations are the same as for 
the £ = 1 case, so the measured phase relations support either 
case equally. Also, since Cm ~ 1/6 for 1 = 2, using the rota- 
tional splitting of the k = 8 and 9 modes yields a splitting of 
2Ckiftat cos f3 ~ 2.2 /iHz, which is close to the measured value 
of 2.06 /iHz. In this case, though, the value of the inclination 
angle, = 22.1°, does not agree with that from the light curve 
fits. Also, the 1 = 2 fits are somewhat less constrained than the 
£ = 1 fits: fitting 500 random amplitudes in the same way as 
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the data, 12% of the fits were better than the fit to the data, 
making the £ = 2 fit significant at the 88% level. 

Even so, given the ability of the oblique pulsator model to 
describe (1) the power in the k=12 region, (2) the phase rela- 
tions between the geometric triplets, and (3) the relative am- 
plitudes within the triplets (to some extent), we believe that 
the above-mentioned inconsistencies do not warrant rejection 
of the model but rather point us in the direction of future im- 
provements. For instance, the amplitude calculations include 
magnetic effects in a fairly simple way; more sophisticated 
treatments may be necessary to adequately model the obser- 
vations. Also, differential rotation in the radial direction may 
be able to resolve the seeming discrepancy of the Cu,i values 
for the £ = 1 case. 

7. DISCUSSION 

In many respects GD 358 is a "simple" star to analyze. 
For instance, the spacing between the periods of its multi- 
plets suggests that we are looking at successive radial over- 
tone numbers of I = 1 modes. This is based both on the fact 
that (1) given its distance, mass, a nd luminosity, only £ = 1 
modes allow a consistent solution (Bra dley & Wingeall994l) 
and (2) in many case s these multiplets are well defined triplets 
(Wing etet al.lll994T) . Furthermore, as a DBV it has no hydro- 
gen layer, so fewer parameters are needed to model its struc- 
ture. Indeed, it is the best studied white dwarf variable and 
signifi cant constraints have been placed on its interior struc- 
ture dMetcalfe. Nather. & Winge3l2^0^lM^tcaife1l200l . 

On the other hand, GD 358 is in many respects a com- 
plicated star to model. For instance, while the amplitudes 
of its modes can be fairly constant during a WET run, over 
timescales of a few months the amplitudes can change sig- 
nificantly. The most extreme case of amplitude change is the 
previously discussed sforzando, in which dramatic amplitude 
changes occurred on a timescale of a day or less. We currently 
have no theory that adequately describes these changes. 

What we do know is that two independent lines of evi- 
dence indicate that GD 358 was hotter during the sforzando 
event. The first is that the average flux of GD 358 relative 
to the comparison stars is larger, and the second is that the 
nonlinear light curve fits indicate a thinner convection zone — 
presumably the result of an increase in surface temperature. 
If we assume that GD 358 remained hotter for approximately 
a day then such a temperature increase would require a total 
energy input of approximately 6 x 10 36 ergs. While detailed 
nonadiabatic models of GD 358's pulsational state just prior 
to the sforzando would be required, rough estimates indicate 
that this amount of energy may typically be present in the 
higher k modes. Thus, if these modes were somehow damped 
and deposited their pulsation energy in the surface layers of 
the star over the period of a day, this would explain GD 358's 
temporary temperature increase. Such a scenario would ex- 
plain both the temperature increase and the disappearance of 
the high k modes, although more detailed models would be 
necessary to check quantitatively the energetics. 

We note that an increase in the apparent amplitude of the 
k = 8 mode is a generic feature of a thinner convection zone 
(e.g., smaller value of to), although the predicted factor (~ 8) 
is still less than what was observed (~ 30). While specula- 
tive, it is also possible that some mechanism may have trans- 
ferred power from the high k to the low k modes, although we 
currently do not understand how this would proceed. Perhaps 
least likely is the possibility that the hotter "equilibrium" state 
of the sforzando allowed mode growth and damping to occur 
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on timescales given by linear theory and that these timescales 
were many orders of magnitude shorter than expected. 

Another possibility is that the event was tied to a magnetic 
phenomenon and that this phenomenon somehow led to a tem- 
perature increase as well as a change in the amplitudes of 
the various modes. There are several reasons that this is not 
completely far-fetched. First, the oblique pulsator model used 
here requires a mechanism to align modes with an axis other 
than the rotation axis, and a magnetic field can provide this. 
Second, the k = 8 and 9 triplets are asymmetric in the sense 
that the m = mode is closer to either the m = +1 or m = -1 
mode (see Figures 18 and 19 of iProvencal etal]|2009l) . Per- 
haps not coincidentally, the asymmetry seen in the splittings 
of the k = 8 and 9 modes reversed themselves at the time of 
the sforzando, and this change has persisted ever since. This 
timescale of years is consistent with wh at is seen in the mag- 
netic cycle of the Sun and other stars (Els worth et al.l[T990l 
Libbrecht & Woodard 1990). Finally, small but definite shifts 
in frequency of order 0.5 /iHz have been seen on timescales 
as short as 3 weeks. These could be a sign of magnetic activ- 
ity: changes in the surface magnetic fie ld could produce slight 
perturbations in the mode frequencies dJones et al.ll 19891) . and 
the timescale of a few weeks for these to take place again 
seems plausible. 

A final unresolved issue for GD 358 is differential rota- 
tion. Taken at face value, the difference in triplet splitting 
between high k and low k modes implies significant differ- 
ential rotation ( Win get et al.lll994j) . In addition, the oblique 
pulsator model as applied in this paper implies differential ro- 
tation: the k = 12 region requires a surface rotation rate of 
~ 2.17 days whereas the k = 8 and 9 splittings give a rota- 
tion rate of ~ 1.5 days. The k = 8 and 9 splittings reflect a 
bulk average of the rotation rate whereas the frequency differ- 
ences within the geometric triplets in the k = 12 region give 
us the rate at the surface. In addition, the frequency split- 
ting between the intrinsic modes in the k = 12 region also im- 
plies differential rotation, albeit with the interior rotating less 
rapidly than the surface. Previously, detailed examinations of 
the frequency splitti ngs in GD 358 as a function of k have 
been inconclusive (Kawaler, Sekii, & Gough l 19991) . but such 
analyses did not take into account the possibility of oblique 
pulsation. 

An intriguing possibility is that the pulsation s themselves 
lead to differential rotation. Townsendj (120091) has recently 
shown that g-modes in massive stars can transport angu- 
lar momentum relatively rapidly compared to evolutionary 
timescales. This transport occurs predominantly in regions 
in which the mode is driven and/or damped and leads to dif- 
ferential rotation. If this effect occurs in white dwarfs then 
it could alter the rotation rate from the surface down past the 
base of the convection zone into the radiative damping lay- 
ers. If the rotation profile of these outer layers is continually 
changing then this would also explain the small shifts in fre- 
quency which have been detected in GD 358. In addition, a 
more dramatic shift in the rotation profile could be associated 
with the sforzando event. Perhaps the outer part of GD 358's 
rotation profile experienced a shift during this event and it has 
persisted in its new state for the past several years. This would 
explain the shift in asymmetry of the k = 8 and 9 t riplets that 
also occurred at this time (see IProvencal et aT1l2009l) . 



8. CONCLUSIONS 

In this paper we have extended our nonlinear light curve fit- 
ting technique to the multiperiodic pulsator GD 358. Our fit to 
the 2006 WET data provides a good match to the light curves 
and we find that the thermal response time of its convection 
zone is tq = 572.9±6.1 s. This is considerably l arger than that 
of the star PG 135 1+489, for which r — 100 s (Mont gomery! 
120051) . This difference in to is consistent with the effec- 
tive temperatures of these stars: the pure He solution for 
PG 135 1+489 yields a T eft which is - 2,000 K hotter than that 
for GD 358. 

We also obtained a fit to the light curve of GD 358 during 
the sforzando event in 1996. These fits showed that GD 358 
had to ~ 42 ± 2 s, a value much less than that determined from 
the 2006 data. This suggests that its effective temperature was 
approximately 2,000 K hotter in 1996 than in 2006 , and th is is 
consistent with the estimate of Wei dner & Ko ester (2003|) that 
the light curve shape suggests T e ff ~ 27 , 000 K. Independent 
evidence of GD 358's brightness relative to comparison stars 
is also consistent with such a temperatu re increase at the time 
of the sforzando (Provencal et al. 2009). The physical origin 
of this temperature increase will be the subject of future work. 

As expected, these data indicate an increase in the 
depth/mass of the convection zone with decreasing T e ff. 
A similar trend is giv en by ML2/a = 1.1 convection 
(Bohm & Cassinelli 1971), although the slope of the theoret- 
ical relation appears less steep than that of the data. In addi- 
tion, lower values of a < 0.6 are excluded. In general, to is 
also a function of logg, albeit a somewhat weaker one. Our 
ultimate goal is to map To as a function of both T e g and logg 
for both the DBV and DAV instability strips. These data will 
provide insight into the physics of convection, still one of the 
largest uncertainties in stellar modeling. They will also serve 
as important constraints for new hydrodynamic simulations of 
convection which are starting to come online. 

Multiple lines of evidence point to some of GD 358's modes 
undergoing oblique pulsation, in particular, the peaks in the 
k = 12 region. First, these peaks can be fit with two sets of 
exactly evenly spaced triplets. Second, the relative phases of 
each of the components within the triplet indicate that each 
originates as a single m = — 1 or +1 mode aligned with the 
magnetic axis; as the star rotates, the magnetic axis precesses 
around the rotation axis, generating a triplet for each intrin- 
sic mode. Finally, the oblique pulsator model qualitatively 
and quantitatively fits the amplitudes of the peaks seen in the 
Fourier transform. Taken together, this marks the first time 
that oblique pulsation has been seen in a white dwarf variable. 

Having now identified the characteristics of oblique pulsa- 
tion in GD 358 we now know what to look for in other white 
dwarf variables; we have found preliminary indications of it 
in other stars and in other data sets of GD 358. As discussed 
in the previous sections, oblique pulsation may prove to be a 
diagnostic of both the magnetic field and its changes as well 
as a diagnostic of differential rotation. This opens an exciting 
chapter in the seismology of these objects. 
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APPENDIX 

EQUATIONS OF OBLIQUE PULSATION FOR £ = 1 AND 2 MODES 

The formulae given below are taken from lUnno et al.l (119891) . Following [Kurtz &" Shibahashi ( 1986) we use u> instead of a for 
the angular frequencies of the modes. We use u/ 0) for the unperturbed frequency of the mode in a non-rotating, non-magnetic star, 
and cjL| for the perturbation to this frequency due to the magnetic field (but not due to rotation). The pulsation axis is assumed 
to be aligned with the magnetic axis, which makes an angle of (3 with the rotation axis, and 9, is taken to be the angle between 
the rotation axis and our line-of-sight. As the pulsation and magnetic axes rotate around the star, we find that a mode of given 
frequency and {£,m} values is split into a triplet of peaks. We tabulate below the resulting time dependence of modes having the 
given values of I and m. 

In the following, x\ is defined to be 
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where C is the rotational splitting coefficient due to the Coriolis force for I = 1 modes and 57 is the angular velocity of rotation at 
the stellar surface. Since x\ is the ratio of the rotational splitting to the splitting induced by the magnetic field it provides a useful 
estimate of the relative importance of the two effects. Also, since it depends on the difference in magnetic splitting between 
m = and \m\ = 1 modes, it is sensitive to the geometry of the magnetic field. With these definitions, the luminosity perturbations 
associated with £ = 1 oblique pulsation are given below: 
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More generally for arbitrary I and m, we define X| m | 
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where C is the rotational splitting coefficient for the appropriate t values. With these definitions, the luminosity perturbations 
associated with 1 = 2 oblique pulsation are given below: 
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